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Recent modern space missions deliver invaluable information about origin of our 
universe, physical processes in the vicinity of black holes and other exotic astro- 
physical objects, stellar dynamics of our galaxy, etc. On the other hand, space 
astrometric missions make it possible to determine with unparalleled precision 
distances to stars and cosmological objects as well as their physical characteristics 
and positions on the celestial sphere. Permanently growing accuracy of space astro- 
nomical observations and the urgent need for adequate data processing algorithms 
require corresponding development of an adequate theory of reference frames along 
with unambiguous description of propagation of light rays from a source of light 
to observer. Such a theory must be based on the Einstein's general relativity and 
account for numerous relativistic effects both in the solar system and outside of 
its boundary. The main features of the relativistic theory of reference frames are 
presented in this work. A hierarchy of the frames is described starting from the 
perturbed cosmological Friedmann-Robertson- Walker metric and going to the ob- 
server's frame through the intermediate barycentric and geocentric frames in the 
solar system. Microarcsecond astrometry and effects of propagation of light rays 
in time-dependent gravitational fields are discussed as well. 



1 Introduction 

The role of high quality reference frames in astronomy has been recognized early by 
both theorists and observers. Astrometric and navigation data rely on observations 
referred to a frame, either local or global. The choice of the reference frame may be 
driven by instrumental considerations or be based upon deeper theoretical grounds. 
This paper deals with the latter subject. 

The original approaches to construct reference frames in astronomy were com- 
pletely based on the concepts of Newtonian gravity and Euclidean absolute space 
and time a. Modern astrometry, however, is operating at the angular resolution 
already exceeding 1 milliarcsecond (see, for example, □) . At this level the primary 
gravitational theory must be the General Relativity Theory (GRT) with a cor- 
responding replacement of the Euclidean space and time by the four-dimensional 
Riemannian space-time manifold. In other words, the theoretical basis of modern 
astrometry must be entirely relativistic. Recognition of this fact is rapidly spread- 
ing in the astrometric community especially after the successful completion of the 
HIPPARCOS missioncl, materialization of the International Coordinate Reference 
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Frame in the sky, development of technologically new ideas in space astrometry, 
and the adoption by the XXIVth General Assembly (GA) of the International As- 
tronomical Union (IAU) of a new relativistic framework for the reference frames 
in the Solar system u. 

In what follows we shall describe in more detail the initial theoretical moti- 
vations that stimulated the development of modern relativistic framework for the 
reference frames in the solar system and new ideas driving the development of the 
theory beyond 2000. Then, in section 2, we give a description of the IAU-1991 
reference frame framework and compare it in section 3 with the present day for- 
mulation adopted by the XXIVth GA of the IAU-2000. A possible way of matching 
the IAU-2000 framework with cosmological model of the universe is considered in 
section 4. The basic theoretical ideas for approaching microarcsecond accuracy in 
relativistic space astrometry is outlined in section 5. 



1.1 Initial Motivations 



Transition from the Newtonian concepts towards more profound theoretical rela- 
tivistic approach for the construction reference frames in modern astrometry and 
celestial mechanics can be traced back to the period 1975-1992 and is associated 
with the breakthrough in the solution of two major problems: development of a 
self-consistent framework for derivation of the post-Newtonian equations of transla- 
±ianal and rotational motion of self-gravitating extended bodies in the solar system 
& □ ancLdpyeJopment of the higher-order relathristic celestial mechanics of binary 
pulsars 0' Er El discovered by Hulse and Taylor E9. The main concern of theorists 
working on relativistic celestial mechanics of the solar system bodies was related to 
the problem of unambiguous interpretation of astronomical measurements and sep- 
aration of small coordinate perturbations from the real physical relativistic effects. 
Later on the question regarding the best choice for the gauge conditions imposed 
on the metric tensor arose, for a number of gauges were used in calculations by 
various groups and there were arguments about their advantages and disadvan- 
tages. But, probably the most serious was the problem of construction of the local 
geocentric (and planetocentric) reference frame in the post-Newtonian approxima- 
tion^). Apparently there existed principal difficulties in solving these problems. 
This was because the simplest approximation of the Earth's gravitational field by 
the Schwarzschild solution could not be considered as accurate enough due to the 
noticeable contribution to the metric from the Earth's rotation and oblateness as 
well as the existence of tidal forces from the Moon, Sun, and other planets. It was 
also recognizedd that the well-known procedure of construction of the Fermi normal 
coordinatestil can not be applied due to the ambiguity in choosing the background 
space-time manifold and deviation of the Earth's center-of-mass world line from 
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geodesic motion. Lack of a rational approach to the problem led to the fact that a 
simple (Euclidean-like) spatial translation of the origin of the barycentric reference 
frame to the Earth's geocenter was used in order to construct a geocentric reference 
frame at the post-Newtonain approximation. However, though such a procedure 
is allowed in GRT because of the coordinate freedom, the Euclidean translation 
to the geocenter produces large relativistic effects having pure coordinate origin 
and, hence, unobservable. In addition, in such a frame the geometric shape of 
the Earth moving along the elliptic orbit undergoes the Lorentz and gravitational 
contractions which must be compensated by spurious internal stresses in order to 
prevent the appearance of unphysical deformations. Similar problems wece also 
met in developing general relativistic celestial mechanics of binary pulsars£3. 

1.2 Motivations in the New Millennium 

Recent technological developments make it necessary to extend the domain of ap- 
plicability of the relativistic theory of reference frames outside the boundaries^ 
the solar system. New generation of astrometric satellites which include-,FAMEt3, 
SEVlli-3, and GAIA (recently adopted as a cornerstone mission of ESAt3) requires 
an absolutely new approach for unambiguous interpretation of astrometric data 
obtained from the on-board optical interferometers. FAME's resolution is about 
100 microarcseconds for stars having 10-th stellar magnitude. The resolution of 
GAIA for the same kind of stars is expected to be already 100 times better. At this 
level of accuracy the problem of propagation of light rays must be basedjXm the 
recently developed post-Minkowskian " Lorentz-covariant" approximation E£l that 
allows us to integrate the equations of light propagation without artificial assump- 
tions about the motion of light-ray-deflecting bodies. Besides, the treatment of the 
parallax, aberration, and proper motion of celestial objects becomes much more 
involved requiring better theoretical definitions of reference frames on a curved 
space-time manifold. 

Additional motivation for improving relativistic theory of reference frames is 
related to the problem of calculation of the incoming signals of the space intcr- 
ferometric gravitational wave detectors like LISA Ell which will consist of three 
satellites flying in space separated by distances of order 1.5 x 10 6 km. Detection 
of gravitational waves can be done only under the condition that all coordinate- 
dependent phenomena are completely understood and subtracted from the signal. 
On the other hand, LISA will fly in the near-zone of the Sun which is considered 
as a source of low-frequency gravitational waves produced by oscillations of its 
interior (known as g-modes)!!^. Motion of the satellites carrying lasers and mirrors 
as well as propagation of light rays along the baseline of LISA in the field of such 
g-modes must be carefully studied in order to provide a correct interpretation of 
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observations. 

Many other relativistic experiments for testing GRT and alternative theories 
of gravity_ako Jpmand advanced theory of relativistic reference frames (see, for 
example, B EM). 

1.3 Historical Remarks 

Historically it was de Sitter B@ who worked out a relativistic approach to build 
reference frames in GRT. He succeeded in the derivation of relativistic equations of 
motion of the solar system bodies and discovered the main post-Newtonian effects 
including gravitomagnetic perihelion precession of a planetary orbit due to the 
angular momentum of the Sun as well as the geodetic precessiotLi'^jS Sitter- Fokker 
effect") that has been verified with the precision of about l%E-jEi These effects 
are essential in the present-day definition of dynamical or kinematical rotation of 
a reference frame u. Later on Lense and Thirring gave a more general treatment 
of the_gravitomagnetic dragging of a satellite orbiting around a massive rotating 
bodycj but hardly believed that the effect can be measured in practice. 

It was GinzburgEH who first realized that the Lense-Thirring effect can indeed 
be measured using artificial satellites and considered the problem of separation of 
coordinate and physical effects. But only xeccntly the experimental verification of 
the Lense-Thirring effect has come about Ej. 

Ginzburg's paper and the success of the soviet space program motivated Brum- 
bergEa to develop a post-Newtonian Hill-Brown theory of motion of the Moon in the 
solar barycentric coordinate system. BaierleinEJ extended Brumberg's approach 
accounting for the eccentricity of Earth's orbit. 

A different original approach to the problem of motion of the Moon and the 
construction of reference frames in GRT was suggested by Mashhoon and Theiss 
in a series of papers (see, for example, l3 and references therein). Instead of 
making use of the post-Newtonian approximations (PNA) they developed a post- 
Schwarzschild treatment of gravitomagnetic effects in the three-body problem. It 
allowed them to discover that the validity of PNA is restricted in time and the 
geodetic and gravitomagnetic precessions are parts of the more general phenom- 
ena involving a long-term relativistic nutation (" Mashhoon- Theiss effect"). The 
same authors also introduced a precise definition of thcJocal geocentric frame with 
the Earth considered as a massive monopolc particleEr c3. 

A more general approach to the_.problem of construction of reference frames 
in GRT was initiated in a (see also E3) , where the matched asymptotic expansion 
technique and a decomposition of gravitational fields of bodies in the (Newtonian) 
multipoles were employed for this purpose. Our approach has been further used for 
the development of the extended Brumbcrg-Kopeikin (BK) formalism for building 
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relativistic astronomical frames. This formalism relies upon Einstein's equations, 
which are solved for the construction of local geocentric and global barycentric 
reference systems, and matching technique, jyhich is used for setting up relativistic 
transformations between these two systems H'E3. 

T. Damour, M. Soffel and C. Xu (DSX) have extended the Brumberg-Kopeikin 
theory by applying the post-Neactonian definitions of the (" Blanchet-Damour" ) 
gravitational multipole momentsLj. Elements of the BK formalism were introduced 
in resolutions of the GA of the IAU in 1991E3. The complete BK-DSX theory is 
presently accepted by the XXIVth GA of the IAU-2000 as a basic framework for 
setting up relativistic time scales and astronomical reference frames in the solar 
system. 



2 The IAU-1991 Reference Systems Framework 

Official transition of the astronomical community from Newtonian positions to 
relativistic concepts began in 1991 when a few recommendations (resolution A4) 
were adopted by the GA of the IAU. 

In the first recommendation, the metric tensor in space-time coordinates (ct, x) 
centered at the barycenter of an ensemble of masses is recommended to be written 
in the form 

900 = _l+^glif) +O(c -4 )5 (1) 

g 0l = 0(c- 3 ) , (2) 
2E7(i,x)" 



9ij = 8. 



■i j 



1 



0(0, (3) 



where c is the speed of light in vacuum, U is the sum of the gravitational potentials 
of the ensemble of masses, and of a tidal potential generated by bodies external 
to the ensemble, the latter vanishing at the barycenter. This recommendation 
recognizes that space-time cannot be described by a single coordinate system. 
The recommended form of the metric tensor can be used not only to describe the 
barycentric celestial reference system (BCRS) of the whole solar system, but also 
to define the geocentric celestial reference system (GCRS) centered on the center 
of mass of the Earth with a suitable function U, now depending upon geocentric 
coordinates. In analogy to the GCRS a corresponding celestial reference system 
can be constructed for any other body of the Solar system. 

In the second recommendation, the origin and orientation of the spatial co- 
ordinate grids for the solar system (BCRS) and for the Earth (GCRS) are de- 
fined. The third recommendation defines TCB (Barycentric Coordinate Time) 
and TCG (Geocentric Coordinate Time) - the time coordinates of the BCRS and 
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GCRS, respectively. The relationship between TCB and TCG is given by a full 
4-dimensional transformation 



TCB - TCG = c- 2 



t 



+ U(t,X E ) ) dt + V E r E 



0(c- 4 ), (4) 



where (x E ) = x E (t) and (v E ) are the barycentric coordinate position and velocity 
of the geocenter, r E = x l —x E with x being the barycentric position of the observer, 
and U(t, x E ) is the Newtonian potential of all solar system bodies evaluated at the 
geocenter apart from that of the Earth. 

In August 2000, a new resolution Bl on reference frames and time scales in 
the solar system was adopted by the XXIVth General Assembly of the IAUQ. 
The resolution is based on the first post-Newtonian approximation of GRT and 
completely abandons the Newtonian point of view on space and time fl. 

3 The IAU-2000 Reference Systems Framework 

3. 1 Conventions for the Barycentric Celestial Reference System 
BCRS is defined mathematically in terms of a metric tensor which reads 

g 00 = -l+ — - — + O{c 5 , 5 

c z c 4 

9oi = — oUX + 0{c- 5 ) , (6) 



c 3 



9ii = % (l + ^wj + 0(e- 4 ) . (7) 

Here, the post-Newtonian gravitational potential w generalizes the usual New- 
tonian potential U and w 1 is the vector potential related with gravitomagnetic 
effects. 

This form of the barycentric metric tensor implies that the barycentric spatial 
coordinates x % satisfy the harmonic gauge condition. The main arguments in favor 
of the harmonic gauge are: (1) tremendous work on GRT has been done with the 
harmonic gauge that was found to be a useful and simplifying gauge for all kinds 
(2) in contrast to the standard post-Newtonian (PN) gauge 
the harmonic gauge can be defined to higher PN-orders, and 
in fact for the exact Einstein theory of gravity. 

Assuming space-time to be asymptotically flat (no gravitational fields exist 
at infinity) in the standard harmonic gauge the post-Newtonian field equations of 

b How the resoltuion Brl— »mg prepared can be traced in the materials of the corresponding 
working groups of the IAUEfl'EZl 



of applications, at 
(see, for example, t 



G 



GRT are solved by 



»(*,«)=<?/ dV + * / dV^,x')|x-x'| , (8) 

w%x)=G [ (9) 
7 |x-x'| 

where cr(i,x) = C - 2 (T 00 + T ss ), cr^x) = c" 1 ^, and T» v = T^(t,x l ) are 
the components of the stress-energy tensor in the barycentric coordinate system, 



3. ,2 Conventions for the Geocentric Celestial Reference System 
GCRS is defined in terms of the geocentric metric tensor 

2W 2W 2 

Goo = -l + ^---^ + 0(c- 5 ), (10) 
c z c 1 

Goa = ~Xw a , (11) 

cr 

G ah = S ab (\ + 1^ + 0(c~ 4 ) . (12) 

Here W = W(T, X) is the post-Newtonian gravitational potential in the geocentric 
system and W a (T, X) is the corresponding vector potential. These geocentric 
potentials should be split into two parts: potentials We and W E arising from the 
gravitational action of the Earth and external parts W cxt and W^ xt due to tidal 
and kinematic effects. The external parts are assumed to vanish at the geocenter 
and admit an expansion into positive powers of X. Explicitly, 

W(T, X) - W E (T, X) + W kin (T, X) + W tidal (T, X) , (13) 
W a (T, X) = W E (T, X) + W& n (T, X) + W&utT, X) . (14) 

The Earth's potentials We and W E are defined in the same way as w E and w E but 
with quantities calculated in the GCRS. Wkin and W£ in are kinematic contributions 
that are linear in X° 



WW = Q b X\ W£ in = - c 2 e abc (Q b ~ n b prcc ) X c , (15) 



where Q b characterizes the deviation of the actual worldline of the origin of the 
GCRS fcpm geodesic motion in the external gravitational field (for more details 
seeMhi) 

Q b = d b w cxt (x E ) - a E + 0(c- 2 ) . (16) 
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Here a l E = dv E / dt is the barycentric acceleration of the origin of the GCRS (geo- 
center). The function f2p roc describes the relativistic precession of dynamically 
nonrotating spatial axes with respect to remote objects: 



prcc 



4 £ vh \ -^v E o k w~ 



t (x fJ ) + 2a fc < t (x B )-i^Q^. (17) 



The three terms on the right-hand side of this equation represent the geodetic, 
Lense-Thirring, and Thomas precessions, respectively. One can prove that flf ncl 
is dominated by geodetic precession amounting to ~ 2" per century plus short- 
periodic terms usually called geodetic nutation. One sees that for Q a = Slp rcc the 
vector potential W^ in vanishes. This implies that dynamical equations of motion 
of a test body, e.g., a satellite orbiting around the Earth, do not contain the 
Coriolis and centrifugal terms, i.e., the local geocentric spatial coordinates X a are 
dynamically non-rotating. For practical reasons, however, the use of kinematically 
non-rotating geocentric coordinates defined by Q a — is recommended. 

Potentials T/(/ tldal and VF* ldal are generalizations of the Newtonian tidal po- 
tential. We also note that the local gravitational potentials W E and W% of the 
Earth are related to the barycentric gravitational potentials we and w E by the 
relativistic transformationstx 



3.3 Transformations between the reference systems 

The coordinate transformations between the BCRS and GCRS are written as 
1 



T = t--[A(t)+v E r E ] 



(18) 



X 1 = rl 



B(t) + B\t) r E + B^ (t) r E r J E + C(t, x) + 0( c - 5 ) 



7 > v E v' E r' E + w eKt (x E )r E + r^o^ri, - \a E r\ 



+ 0(c- 4 ), (19) 



where functions A(t), B(t), B l (t), B^ (i),C(i,x) can be found inB. Let us also re- 
mark that the harmonic gauge condition does not fix the function C(t, x) uniquely. 
However, we prefer to fix it in the time transformation for practical reasons. 



4 Matching the IAU-2000 Framework with Cosmological Reference 
Frame 

The rapidly growing accuracy of astronomical measurements makes it necessary 
to take into account some important cosmological effects for an adequate inter- 
pretation of optical and radio observations of cosmological lenses, anisotropy in 
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cosmic microwave background radiation, etc. For this reason the matching of the 
IAU-2000 framework for reference systems in the solar system with the cosmologi- 
cal reference frame becomes vitally important. In this section we outline the main 
ideas in this matching for the spatially flat Friedmann-Robertson- Walker (FRW) 
cosmological model. 

First, we consider the perturbation h a p of the gravitational field as 

g a! i = a 2 (T){7 laP + h a i ) ) , T*fi = -h«P + ± v a f , h, h = V af3 h af3 , (20) 

and impose the quasi-harmonic cosmological gauge conditions [] 

dp-f* = -2H (V° - ^ q °t) , H{t) = - a , (21) 

which eliminate almost all first order derivatives of the metric perturbation. Then, 
the linearized Einstein equations read 

□ 7 «/3 _ 2H d Tl al3 + 2{2H + H 2 )h af3 + A(H - 2H 2 )rf^ a 'h^ (22) 

- 2(H - H 2 ) V af) h ao = \&na 4 8T af3 , 

where a dot over "the Hubble parameter" H denotes the time derivative and ST a @ 
is the tensor of energy- momentum of perturbing source (" solar system" ) . 

For the particular case of matter-dominated background FRW cosmological 
model the reduced linearized Einstein equations (p2|) read 



a gW - 8H 2 w = -47ra 4 (r) <5(T 00 + T kk ) , (23) 

□ ff ( 7 + \w) = -4™ 4 (r) S(T™ - 7 -T) , (24) 

a gW l - 5ffV = -47ra 4 (r) ST n , (25) 

U gW l] = -47ra 4 (r) 5T <ij> , (26) 

where w = - (7 00 + j kk ) /4 , w l = -7 0i /4, w 13 = _ T <^>/4fJ and H(t) = 2/r. 
For the differential operator D g one has 

d 2 4 d 
dr 2 t dr 



so that the solution of the inhomogeneous equation U g F{r, x) = — A-kcl (t)ST(t, x), 
can be found making use of the replacement ELI 

, , 1 d /*(r,x)\ , x 

F(r,x =-5- ■ (28) 



r 



dr 



c These conditions were also independently discovered in the papcr|23 for the case of de Sitter 
space-time. 1—1 
d The square brackets around spatial indices denote symmetric and trace-free tensor 
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transforming the equation in question into 



\l T (^^) = -W(r)*r(T,x) , = A - 81 , (29) 

which has a simple particular solution 

*(t,x)=t J r^^Tj ^ sa 4 {s)ST{s, X r )ds- j d 3 x' ^ aa 4 (s)(ST(s, x')ds , 

^ ' 1 (30) 

where t' = t — |x — x'| is a retarded time in a flat space-time. 

Matching this solution with that defined in the BCRS of the 1AU-2000 frame- 
work as flat-space retarded potentials is achieved after the replacement r — 1 + Hgt 
and expanding all quantities depending on r in the neighbourhood of the present 
epoch t = 1 along with making use of the tensor transformation law for metric 
tensor. Results of the matching will be described elsewhere. 

5 Relativistic Microarcsecond Astrometry 

The IAU-2000 framework for reference frames requires new advanced theory of 
astrometric data analysis. The key issue in this theory is solution of the problem 
of propagation of light rays with account for as many relativistic effects as required 
for making unambigious interpretation of the data. Previous approaches for calcu- 
lating light ray propagation in the framework of relativistic astrometry were based 
on making use either metric tensor of exact solutions of GRT (Schwarzschild, Kerr, 
etc.), or metric tensor of the post-Newtonian approximation (PNA), or metric ten- 
sor of the plane weak quadrupolar gravitational wave. All these approaches have 
difficulties and/or inconsistencies in describing light propagation at the microarc- 
second threshold. The question arises how to avoid these difficulties and what is 
the proper formalism to deal with the problem. Treatment of this problem has 
been given mcexilly in the framework of the " Lorentz-covariant" theory of light 
propagation lj ufj. and the formalism developed has an unrestricted ability to 
make calculations at any desired order of approximation with respect to the small 
parameter v a /c as well as many other advantages making it a powerful tool for 
theoretical predicitons of various relativistie-ieffects being detectable by modern 
astrometric technique (see, for example, E-irEj). 
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